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Abstract 

The finite-temperature effective potential customarily employed to describe the physics of cosmo- 
logical phase transitions often relies on specific gauge choices, and is manifestly not gauge-invariant 
at finite order in its perturbative expansion. As a result, quantities relevant for the calculation of 
the spectrum of stochastic gravity waves resulting from bubble collisions in first-order phase tran- 
sitions are also not gauge-invariant. We assess the quantitative impact of this gauge-dependence 
on key quantities entering predictions for gravity waves from first order cosmological phase transi- 
tions. We resort to a simple abelian Higgs model, and discuss the case of gauges. By comparing 
with results obtained using a gauge-invariant Hamiltonian formalism, we show that the choice of 
gauge can have a dramatic effect on theoretical predictions for the normalization and shape of the 
expected gravity wave spectrum. We also analyze the impact of resumming higher-order contribu- 
tions as needed to maintain the validity of the perturbative expansion, and show that doing so can 
suppress the amplitude of the spectrum by an order of magnitude or more. We comment on open 
issues and possible strategies for carrying out "daisy resummed" gauge invariant computations in 
non-Abelian models for which a gauge-invariant Hamiltonian formalism is not presently available. 
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I. INTRODUCTION 



The search for gravitational waves is entering an exciting phase. The current generation of 
experiments is aheady dehvering interesting results, including recent limits on the amplitude 
of stochastic gravitational wave backgrounds from the LIGO Collaboration [U [2]. Rapid 
advances in the development of space-borne detectors that might be operational in the 
relatively near future are also expected 

As pointed out long ago [5J, cosmological phase transitions in the early universe might 
produce an imprint in the form of a stochastic background of relic gravity waves. These would 
arise as a result of the collision or turbulent motion of bubbles of "true vacuum" expanding 
and eventually filling the metastable vacuum in a cosmological first order phase transition [6] . 
The resulting signal might be large enough to be detectable by the next generation gravity 
wave search experiments, providing a unique window on the early cosmological history of 
the universe 

The spectrum of the gravity wave background arising from a first order cosmological 
phase transition is controlled by two physical properties of the phase transition itself: (a) 
the amount of false vacuum energy liberated during the phase transition - in other words, the 
latent heat associated with the transition; and (b) the bubble nucleation time scale, which 
gives a measure of how rapidly the phase transition occurs relative to the early universe 
Hubble expansion rate [HI E] • The resulting gravity wave normalization and spectral peak 
can be estimated as a function of these two physical quantities. Detailed analytical [Hlin] and 
numerical [101 E] studies exist that relate the two parameters to the predicted spectrum, in 
particular for the case of detonations, where the speed of the bubble wall is larger than the 
speed of sound (see Ref. [11] for a discussion of the opposite case of deflagration). 

One class of models where a strongly first order phase transition is a necessary ingredient 
is electroweak baryogenesis |[12j . In the presence of B- violating electroweak sphalerons in the 
Standard Model(SM) and sources of CP-violation beyond those of the CKM CP violating 
phase, the electroweak phase transition (EWPT) can produce a sufficiently large baryon 
number density to explain the observed baryon asymmetry (for recent studies see e.g. [T3V - 
[25] . To prevent the washout of the produced asymmetry, the phase transition must be 
strongly first order, thus necessarily producing gravity waves. Interestingly, the typical 
frequencies at which gravity waves from the EWPT are red-shifted today falls in the milli- 
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Hertz to Hertz range, where the expected sensitivity of the space-based interferometer LISA 
is maximal. The search for a gravity wave rehc from the EWPT is therefore especially 
intriguing and promising (see e.g. f2E[ f27\ ). 

From a particle physics perspective, the determination of the details of an EWPT de- 
pends on the calculation of the finite-temperature effective action Fefr as a function of the 
background bosonic fields — denoted generically here as ip{x) — that are present in the the- 
ory. In the case of the SM, lattice computations and the LEP lower bound on the mass of 
the Higgs boson imply that electroweak symmetry breaking in a SM universe occurs through 
a cross-over transition |28t |2S]. To obtain a strongly first-order EWPT as needed for both 
electroweak baryogenesis (EWB) and the associated relic gravity waves, one must augment 
the SM scalar sector by the addition of new scalar fields, such as a second Higgs SU(2)i: 
doublet as in the Minimal Supersymmetric Standard Model (MSSM) p^HlQ] or a real singlet 
in minimal extensions of the SM scalar sector (see e.g. [2Q]). Electroweak bubble nucleation 
occurs when a combination of the one of these new fields and the neutral component of 
the SM Higgs doublet becomes non-zero. The properties of relic gravity waves produced by 
bubble collisions then follows from the behavior of rcff[<y5(x)]. 

The most theoretically robust computations of the effective action are performed using 
non-pertrubative (lattice) methods. Given the cost of such computations, however, this ap- 
proach is not feasible for exploring EWPT dynamics in a variety of beyond the Standard 
Model (BSM) scenarios. Consequently, one must resort to perturbation theory which, in 
turn, requires introduction of gauge-fixing. As recognized long ago [30], perturbative com- 
putations of the effective potential — and more generally refr[v?(a;)] — generically lead to a 
gauge- dependent function. Physical quantities like the latent heat or the bubble nucleation 
rate should not, however, exhibit any gauge dependence. 

In fact, general arguments imply that the critical temperature [3Tj and the bubble nucle- 
ation rate [32] are gauge independent. These statements follow from the so-called Nielsen 
identities and their generalization [33l|3l] that describe the dependence of the effective action 
on the gauge fixing condition imposed on the quantized fields. In particular, the effective 
action refr[(y5(a;)] is gauge-invariant when (p{x) is an extremal configuration, that is, one sat- 
isfying the equations of motion^. In principle, then, one should be able to obtain physical, 

^ In the case of the effective potential, (p{x) = const = ipnun is just the spacetinie independent background 
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gauge-invariant quantities relevant to the EWPT from Tcs by working with an appropriate 
set of extremal field configurations. 

In practice, a non-trivial gauge-dependence can enter perturbative computations from 
an inconsistent truncation of the perturbative expansion [35j. In the context of sphaleron- 
induced baryon number washout, following on earlier work by Refs. [36H38] . it was shown in 
Ref. [35] how a consistent, systematic order-by-order approach can yield a gauge-invariant 
perturbative result. However, most of the remaining literature on the topic of the EWPT, in- 
cluding the context of gravity wave production and of baryogenesis relevant here, appears to 
suffer from gauge-dependence (typically, perturbative calculations of EWPT-relevant quan- 
tities have been performed in the Landau gauge)^. Apart from the point of principle, the 
question then arises as to the quantitative impact that this gauge artifact has upon predic- 
tions of observable quantities. 

In the present study, we address this question as it pertains to computations of gravita- 
tional wave spectra from a first order EWPT. To that end, we consider the simplest model 
involving scalar fields charged under a gauge group: the Abelian Higgs model, also known as 
the Coleman- Weinberg or scalar QED model. We then resort to a class of gauges known as 

(or renormalizable) gauges, and we calculate the effective potential at finite temperature 
including its explicit dependence on the parameter ^. 

Studying the Abelian Higgs model has two clear advantages. First, its parameter space 
is small and easily analyzed. Second, and more importantly, one may calculate its effective 
potential using a gauge- invariant Hamiltonian approach [32] whose results can be compared 
with those obtained from the gauge-dependent approach. We then calculate quantities rele- 
vant to the character, strength and duration of the EWPT both in the gauges (including 
Landau gauge) and in the gauge-invariant approach, and systematically compare the results. 

We find that the gauge choice may have a dramatic impact, amounting to several orders 
of magnitude, on the inferred gravity wave spectrum, and even on the first- or second-order 
character of the phase transition itself. We also observe that the Landau gauge results closely 
matches, at least for the Abelian Higgs model, the results using the explicitly gauge-invariant 
Hamiltonian formulation. This situation is perhaps not surprising, given the arguments in 

field {e.g., Higgs vacuum expectation value) that gives a minimum or maximum of the potential. 
^ To our knowledge, there exist no non-perturbative computations of all of the quantities relevant for 
predictions of GW spectra. 
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Ref. [10] (see below). We caution however that this conclusion might not be easily gener- 
alizable to non-Abelian gauge theories and that even in the Abelian Higgs model, a gauge- 
invariant resummation of higher-order terms that would otherwise spoil the convergence of 
the perturbative expansion remains to be developed^. To underscore the importance of the 
latter problem, we study the impact of including the "daisy resummation" in gauges. We 
find that in the Landau gauge, inclusion of the resummation typically reduces the overall 
amplitude of the GW spectrum compared to the gauge-invariant but un-resummed result. 
We then comment on strategies to tackle these issues in non-Abelian models (such as the 
Standard Model or its supersymmetric extensions) - including those of Ref. [35] • 

We begin in section |TT] with the explicit calculation of the gauge-dependence of the effective 
potential, at both zero and finite temperature, for the Abelian Higgs model. In section \LU\ 
we explain the calculations required to predict gravitational wave spectra and other physical 
observables related to the phase transition. Finally, we present our results and conclusions 
in sections [IV] and |Vl 

II. GAUGE DEPENDENCE OF THE EFFECTIVE POTENTIAL 

We are concerned here with an Abelian Higgs model encompassing a complex singlet 
Higgs field with Lagrangian 



where = df^Ai, — d^jA^ and = d^ — ieA^ are the standard electromagnetic tensor and 
covariant derivative, respectively. The potential Vq is 



The tree- level vacuum expectation value is = and the bare Higgs mass at the vacuum 
expectation value (vev) is m\ = 2m^. 

In order to perform perturbative calculations, we must add gauge-fixing and ghost terms 
to the Lagrangian. In the gauge, these are 




(1) 




(2) 



1 



(3) 



■gf ~l~ ^ghost 



2^ 



A gauge-invariant prescription for estimating these terms was developed in Ref. |35j . 
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where c is the Grassmann-valued ghost field, and where we have spht into its real and 
imaginary components: (p = a + irj. We choose the vev such that (r/) = 0, making a the 
physical Higgs boson and rj the non-physical goldstone boson. Note that in Landau gauge 

= 0, fixing dfj^A'^ = 0), the ghost field completely decouples from the theory. 

To include finite-temperature corrections, we must go to (at least) one-loop order in the 
efi^ective potential. At zero-temperature, the calculation of the one-loop effective poten- 
tial proceeds by taking the trace of the inverse propagators for each particle. This yields 
terms like ^ / log(A;^ + mj(a)), although determining the proper expression for the ^- 
dependence of the gauge boson is somewhat complicated by the need to sum over Lorentz 
indices ( see, e.g., Ref. [S]). Using MS renormalization, the full one-loop zero-temperature 
potential is 



where rii are the degrees of freedom for each particle and c = 1 for the gauge boson's 



particle masses and their degrees of freedom. Several of the masses are gauge-dependent, 
and, precisely because of this fact, the effective potential is also gauge-dependent. Note 
that at both the origin and the tree-level vev (a^ = = 2m?' /X) the gauge-dependence 
disappears [35], as expected from the Nielsen identities [331 [M]- However, the value of a that 
minimizes the one- loop effective potential is not gauge invariant. 

For the particular case of the Abelian Higgs model, Fischler and Brout [30] defined 
an effective potential from the vacuum-to-vacuum ^-matrix element without resorting to 
the introduction of sources, a procedure that contrasts with the conventional definition in 
terms of the Legendre transform of the source-dependent generating functional, Z[j]. In 
this context, the "free-energy" is minimized by a spacetime-independent background field 
only in the Landau gauge, whereas in other gauges the minimizing fields must carry a 
spacetime dependence. Consequently, only in the Landau gauge does the minimum of the 
effective potential in the absence of sources characterize the presence or absence of symmetry- 
breaking. For the formulation with sources, a spacetime-independent background field will 

^ In the literature, authors generally use c = | for all three physical modes of the gauge boson. This only 
makes a difference if one includes thermal masses in the zero-temperature potential, which is a small 
correction that few authors include. 





transverse modes and c = | for its other modes and all other particles^. Table □ lists all 
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particles 


d.o.f. 


(masses)^ 


(thermal masses)^ 


transverse gauge polarization 


2 






longitudinal gauge polarization 


1 






time-like gauge polarization 


1 






higgs boson 


1 




(iA + ie2)T2 


goldstone boson 


1 


-m2 + (iA + ee2)c72 




ghosts 


-2 


eeV2 





TABLE I: Particle content of the Abelian Higgs model, including Fadeev-Popov ghosts. One 
ghost effectively cancels the contribution from the unphysical time-like polarization, while the 
other cancels either the longitudinal polarization (at o" = 0) or the goldstone boson (at a = v). 

yield the minimum of energy in any gauge of the form given in Eq. ([s]), implying equal values 
of the minima of the effective potential for any choice of ^ [33l |3l]. While the formulation of 
Ref. jlQ] is manifestly gauge invariant by construction, its relationship with the development 
in terms of sources has not to our knowledge been clarified. That being said, the arguments 
of Ref. [lO] are suggestive that results obtained with the Landau gauge effective potential 
may be most physically reasonable. Indeed, we find close numerical agreement between 
Landau gauge quantities and those obtained using an explicitly gauge-invariant Hamiltonian 
formalism (see below). We emphasize, however, that this agreement does not carry over to 
the non- Abelian case. 

The finite temperature contribution can be derived similarly to the zero-temperature 
contribution, except that the integral over momenta is replaced with a sum over Matsubara 
modes: / dk^ ^ J^p -)■ This yields 

Vi(a,r>o) = ^ E J 

particles 

where 

J(x') = dy log (\ - e- . (6) 
In the high-temperature (low-x) limit, 

J(x2) ^ + - -a;3 - - log - - 0(x6) (7) 
^ ^ 45 12 6 32 ^ Ob ^ ^ ^ ' 

where loga^ = | — 27^; + 21og(47r) and is the Euler constant [12]. All higher-order terms 
are simple polynomials in . Again, the gauge dependence disappears at = and at 
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(5) 



(T = 0, but it is non-trivial everywhere else. We also observe that the general arguments in 
Ref. |im do not depend on whether one works with a Minkowski or Euclidean formulation 
of the functional integral appearing in the generating functional, so that even at finite-T use 
of the Landau gauge is equivalent to a gauge-invariant formulation for the Abelian Higgs 
model. 

A. Thermal Mass Corrections 

It is well-known that near the critical temperature for a phase transition, validity of the 
perturbative expansion of the effective potential breaks down. Quadratically divergent con- 
tributions from non-zero Matsubara modes must be re-summed through inclusion of thermal 
masses in the one-loop propagators [131111]: m^(cr) — t- mljj^a) = {a) + m^j^^^^{T) . Tablejl] 
lists all thermal mass corrections (see ref. |1S] for further discussion and explicit calculations 
of the masses). Generally, one performs this "daisy resummation" by only including the 
thermal masses in the zero-mode propagators, which results in a mass correction to only the 
cubic term in the effective potential. It is slightly more convenient from a computational 
standpoint to include the corrections in all propagators, although we do check that this only 
makes a small difference in the resulting potential. 

B. Alternative Gauge Invariant Formulation 

For comparison, we also examine the gauge-invariant effective potential put forward by 
Boyanovsky et al. [3S]. These authors derive the potential by working in the Hamiltonian 
formalism using only gauge-invariant physical states. In this case, there exist only four 
independent degrees of freedom (two transverse gauge, one longitudinal gauge, and Higgs), 
with no need for ghost cancellations. The unrenormalized one-loop effective potential is 



where uJ^ = k'^ + m\, and = A;^ + m\ arise from the transverse gauge and Higgs degrees of 



contribution of both the gauge boson's longitudinal polarization and the Goldstone boson. 
The order parameter x is a spacetime- independent, gauge- independent shift of the field, and 




(8) 



freedom, respectively, and the plasma frequency 



(/c^ + m^g){k'^ + m'^)/k'^ contains the 
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the gauge, Higgs, and Goldstone masses m^, m^, and rUg are given in table I with ^ = 
and 0" — )■ X- The tree-level potential is unchanged. 

The first two contributions to Vi{x) exactly match the transverse gauge polarization 
and Higgs boson contributions to the potential in gauge. At the tree-level extrema, 
the plasma frequency also matches the contributions from all other modes. However, away 
from the tree-level extrema, the plasma frequency does not match and looks similar only 
to Landau gauge = 0), as one would expect from the general arguments in Ref. |10]. 
Therefore, we anticipate the Landau gauge will provide a close approximation to the gauge- 
independent Hamiltonian result, a conclusion similar to what was found in Boyanovsky et 
al. 1391. 



Using MS regularization (see the appendix), we find that the plasma frequency contri- 
bution to the one-loop potential is 



V,p{x,T>0) 



647r2 

rp4 



{m% - rr?^f ( log 



2 2 

— m 



+ Arrirpm. 



27r2 



2^2 ' 2^2 



where J2 is calculated by Boyanovsky et al. to be 



dy y"^ log 



1 



(9) 
(10) 

(11) 



In the high temperature expansion, Boyanovsky et al. find that their gauge invariant poten- 
tial is the same as the potential in Landau gauge up to the cubic terms, but the equality 
breaks down beyond this. 

In what follows, we compare results for the GW spectra using the Lagrangian and gauge- 
invariant Hamiltonian methods. We observe that the daisy resummation of higher-order 
contributions was not considered in Ref. [39] , and it is not immediately clear how one would 
do so. Consequently, when comparing results in the two approaches, we will not include the 
Daisy resummation in the Lagrangian formulation. 



III. CALCULATION OF GW PARAMETERS 



We calculate several parameters of interest to gravitational wave production from early 
universe phase transitions using the Abelian Higgs model with full and explicit gauge depen- 
dence in the class of gauges. These include the transition temperature T^,, the minima of 
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the low and high-temperature phases at the transition, the relative change in energy density 
a, and the approximate duration of the phase transition 



A. Calculating the Transition Temperature 

In first-order cosmological phase transitions, the low-temperature phase develops by nu- 
cleating bubbles within the high-temperature phase (see Refs. |l6l ST] for original work on 
cosmological transitions). A critical bubble — one whose surface tension exactly balances its 
outward pressure — is given by the 0{3) symmetric action 



S-^ = An { r dr 



oo 

2, 



V^V + mr),T) 



(12) 



2\dr J 

subject to the constraints that the field is smooth at r = and in the high-temperature 

minimum at r = oo. Smaller bubbles collapse, while larger bubbles grow and eventually fill 

the universe with the new phase. Equation [12] yields the radial equation of motion 

d'^a 2 da d , , 

dr"^ r dr da ' ' 

which we solve using the 'undershoot /overshoot' method (see e.g. Ref. [H]). 

To find the exact transition temperature T^,, we must determine when the low-temperature 
phase nucleates at least one bubble per Hubble volume. The nucleation rate goes roughly 
as r oc T^e""^^/^, where the constant of proportionality can be found largely on dimensional 
grounds. For electroweak scales, this gives a transition temperature determined by S'3/T* ~ 
140 (see e.g. Ref. jUj). Note that the exponent changes very rapidly, so determining the 
exact form of the coefficient is quantitatively unimportant. 

Finding the minima of the low- and high-temperature phases can be a non-trivial task, 
especially since the high-temperature minimum is not necessarily at cr = and intermediate 



minima can develop for ^ > (see section IV below). Our strategy, however, is fairly 



straightforward. We first observe that the transition occurs in the range Tc > > T^in, 
where Tmin is the lowest temperature at which the original, high-temperature phase exists 
and Tq is the temperature at which the minima of the potential in the two phases are 
degenerate. We then trace the low-temperature minimum upwards from T = and the 
high-temperature minimum downwards from T = Tq by numerically integrating 

dcTrn.in 



dT 
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At each point, we calculate the transition rate by finding 5*3. Following the evolution of the 
minima and Sc we then determine the temperature at which S3/T = 140. 



B. Calculating the Latent Heat and Transition Duration 

With the transition temperature in hand, the relative change in energy densities a and 
transition duration /3~^ follow without much effort. When evaluated at its minimum, the 
effective potential is the same as the free energy density of the system^. Therefore, the 
energy density difference between the two phases is 

Ap = [V{(rhot, %) + ShotT*] - \V{acoid, %) + Scoid%], (15) 

where the entropy density is s = —dV/dT. Note that at = Tc, Ap is identical to the 
latent heat. The quantity of interest in the production of gravitational waves is a = Ap/ prad, 
where prad = ^i^T'^ and g^, is the number of relativistic degrees of freedom at the phase 
transition, which we take to be 100. 

Writing the bubble nucleation rate as F = FqC^*, gives the approximate phase tran- 
sition duration. For a radiation dominated universe, 



/3 y d{Ss/T) 



(16) 



dT 

where is the Hubble expansion rate during the transition (see e.g. Ref. [H]). 



C. Calculating GW Spectra 

We employ here the analytical approximation provided in ref. to the numerical simu- 
lations carried out in that same work. We refer the Reader to ref. [H] for further insights on 
the results of ref. [10] The gravity wave spectrum (more precisely, the gravity wave energy 
density per frequency octave) from collisions at production is parameterized by 

^GW*{f*) = ^GW* (17) 

^ Here, we neglect kinetic energy contributions associated with non-vanishing gradients of the background 
field in the bubble walls. 



12 



where the two exponents, obtained from fits to the numerical results, are set to a = 2.8 and 
b = 1.0. The spectrum is redshifted according to 



/inn\V3 

h^riGw = 1-67 X 10-^ ^Gw* f j (19) 

with the functions and A depending on the bubble wall velocity Vb (and hence implic- 

itly on the relative energy density difference a) according to the following parameterization 
(again as given in ref. [lO]): 

Finally, we employ the following parameterization for the bubble wall velocity |19] 



V, = — , (23) 

valid in the limit of interest here of strongly first order phase transitions. Note that the 
overall amplitude scales as h'^Qcw oc g* ^^'^ for a ^ 1, so it can be changed by several orders 
of magnitude by choosing a model with a different g^,. 



IV. NUMERICAL RESULTS 



In the most basic model without any additional fields, four parameters determine the 
effective potential: the tree-level Higgs mass rrih, the tree-level vev v, the gauge coupling e^, 
and the renormalization scale /i. We vary only the first two of these, keeping v = 246 GeV 
and fi = 1 TeV fixed. We include phase transition calculations using the gauge-invariant 
Hamiltonian formalism of Boyanovsky et al. |39] without resummation, shown in the figures 
as solid arrows. 

In order to generate a fairly strong first-order phase transition, the gauge boson mass must 
be relatively large. However, if the mass is too large then the one-loop zero-temperature 
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potential overwhelms the tree-level potential and perturbation theory is unreliable. At the 
tree-level vev, Vo{a = v) = -|Aw^ and Vi{a = v,T = 0) = ^{e^v^f[\og{^) - |] plus a 
small contribution from the Higgs boson. These two are roughly equal when ~ 4 A or 
g2 ^ 2?^^ slightly more conservative, we demand that < 

Figures |2]-[4] display our results for the gauge dependence of the different phase transition 
properties. All four properties — the transition temperature, the values of a corresponding to 
the minima of the phases, the relative change in energy density, and the transition duration — 
heavily depend upon the choice of gauge. We also show the impact of including the daisy 
resummation, as discussed above. 

Three broad features emerge from these figures. The results obtained with the Hamil- 
tonian formulation most closely match the results of Landau gauge = 0). However, 
the match is not exact. Most significantly, the Hamiltonian approach yields a small but 
measurable increase in the critical and transition temperatures. 

Second, at ^ = 0, there can exist significant shifts in the GW-wave relevant parameters 
due to the inclusion of the daisy resummation. Generally, one finds that the values of a 
are decreased while are increased due to the inclusion of the resummation, implying a 

reduced amplitude and higher peak frequency in the GW spectrum. This significant depar- 
ture from the fully gauge-invariant results (albeit within the Landau gauge) suggests that 
developing a gauge-invariant daisy resummation procedure will be essential for obtaining 
physically realistic predictions. 

Third, the dependence on can both exacerbate these differences and lead to new phase 
structures that are clearly unphysical artifacts of the gauge choice. For example, for rrih = 
120 GeV (fig.|4]), the phase transition becomes second-order at high ^, so the change in energy 
density goes to zero and /3 goes to infinity. A change in gauge can also lead to secondary 
minima and secondary transitions (see fig. [T|. In figs. [2]-[4| we always perform calculations 
for the transition with the largest change in vev, even when this transition happens after 
initial symmetry breaking. This leads to the discontinuities in figs. [3] and |4j Given the 
unphysical nature of these artifacts, we do not discuss them further but simply point out 
the danger in this context of attempting to draw physical inferences from a gauge- dependent 
calculation. 

Finally, we present our calculations for various gravitational wave spectra in figures |5 - 10 



We make comparisons of the Hamiltonian approach and gauges without daisy resumma- 
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FIG. 1: An example of multiple phase transitions in the same model. Here, rrih = 35 GeV, 
= and ^ = 3. Since the existence of the secondary phase transition is gauge-dependent, it 
is clearly non-physical. 



tion in figures [5]-[7[ and include the effects of resummation in figures [8 -10 Again, Landau 
gauge and the Hamiltonian approach produce very similar results. However, a change in the 
gauge parameter produces very large changes in both the calculated amphtude and peak 
frequency of the wave. Without daisy resummation it does not appear possible, or at least 
feasible, to determine which way this change will manifest without doing the full numeri- 
cal calculation. With resummation, an increase in the gauge parameter tends to make the 
phase transition less strongly first-order, thereby decreasing amplitude and increasing the 
peak frequency of the resulting spectrum. 



V. DISCUSSION 



We have thus far presented numerical calculations of strongly first-order phase transitions 
and spectra of the resultant gravitational waves in the Abelian Higgs model both for various 
values of the gauge parameter ^ and in two gauge-invariant formalisms. The gauge-invariant 
Hamiltonian formalism closely matches Landau gauge. We find that small changes in ^ can 
produce large changes in calculated physical quantities, implying that attention to gauge 
invariance in GW computations is essential for reaching physically meaningful predictions. 
Moreover, we find that in the Lagrangian formalism, the result may be significantly affected 
by inclusion of the daisy resummation, a conclusion similar to what has been observed in 
the context of sphaleron rate computations . 
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From these observations, we conclude that the use of a non-gauge-invariant framework and 
the neglect of daisy resummations in computations of GW spectra for non-Abelian phase 
transitions are likely to lead to physically unreliable predictions. At present, it appears 
that the generalization of gauge-invariant perturbative methods applicable in the Abelian 
Higgs model to non-Abelian spontaneously broken theories is not straightforward. The 
Hamiltonian formalism does not easily carry over to the non-Abelian case and, given how 
drastically observables change with a change in gauge parameter, one should not trust gauge- 
dependent calculations for anything other than rough estimates. 

We see several possible directions. First, one can compute thermodynamic quantities 
of interest (such as a and T^) as well as the bounce action with Monte Carlo methods, 
thereby circumventing the gauge problem at the outset while including all higher-order 
effects (including those entering daisy resummed perturbation theory) by construction. The 
results would undoubtedly be the most reliable theoretically, but this approach is unlikely 
to be practically feasible for surveying a wide variety of models or exploring wide regions of 
parameter space for models like the MSSM. As an alternative, following Ref. o^ig can use 
-Rg gauge and the Nielsen identities to ensure gauge-independence at each order in h. The 
latter approach is relatively straightforward conceptually, but computationally involved, 
as one must go to at least second order in the loop expansion. It appears particularly 
challenging in the case of the tunneling rate computation. From a more formal side, it may 
be possible to construct a gauge-invariant Hamiltonian formalism for spontaneously-broken 
non-Abelian gauge theories. Although we are not aware of any work in this particular 
direction, we note that such a formulation has been achieved in the absence of spontaneous 
symmetry breaking for the specific case of quantum chromodynamics (see, e.g.^ Refs. [50|[5T] 
and references therein). 
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Appendix A: Calculating 1-loop zero-T potential 

To calculate the one-loop potential, we must examine integrals like J (f^cu. For cu^ — 
k"^ + m^, this gives the standard one- loop potential associated with a particle of mass m. 
However, in the gauge invariant approach of Boyanovsky et al. the plasma frequency has the 
form u!p — {k^ + a){k'^ + /3)/k'^. They perform the integral using a cutoff regulator, but we 
would like to use dimensional regularization in order to better compare with the gauge. 

The potential associated with the plasma mode is given by 

" 2 J {2t,Y V(A;2 + a)(A;2 + /3)J ^ ^ 

with d — and n — and is a mass dimension that balances the integration measure. 
For n — I > 0, the integral converges. Performing the integral over the d-dimensional sphere 
yields 

K^ .. / dkk'-' ^ 



P (47r)'^/2r(rf/2) y \{k^ + a){k^ + 

^ ^ ^ dk p"/^-^ ( ) . (A2) 



{AnYnV{d/2)2 J ^ \{p + a){p + P) 
We can introduce a Feynman parameter to rewrite the fraction as 



Using this, and the definition of the beta function 

' dx x'^-\l - xf-' = B{a, h) = Y^^, (A4) 
one can show that 

Then, using the generalized binomial theorem, 

^ 1 T{n-d/2) ^ r(d/2 - n + l)r(d/2 - l)T{n + I) g^.^.i i 
" (477)^^/2 2r(d/2)r(n)2 r(d/2-n + i-0r(/ + i) ^' ^ ' 
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where we demand that |/3| < |a|. Expanding this out in e 



3-d 
2 ' 



one finds 



647r2 



1 




+ -fE- log(47r) 




a — /3 3 
~f? 2 



) 



+ 4af3 +0(e), (A7) 



where is the Euler-Mascheroni constant. In MS regularization, we simply subtract out 
the term containing 1/e, as well as the 'Je and log(47r) terms. Note that for /3 = 0, this 
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FIG. 2: Calculated gauge dependence of phase transition parameters for a low-mass Higgs boson. 
In all panels, black (grey) lines denote models with (without) resummation. The arrows denote 
values corresponding to the solid lines, but calculated in the gauge-invariant Hamiltonian formal- 
ism. All quantities along the y-axes are in units of GcV, except for P/H which is unitless. In 
the first panel, solid, dashed and dotted lines denote the transition temperature T*, the critical 
temperature Tc, and the minimum temperature at which the hot phase exists. In the second panel, 
solid and dashed lines denote the minima of the cold and hot phases. The third panel shows the 
relative difference in energy densities at both the critical temperature (dashed line) and the actual 
transition temperature (solid line). The final panel gives (3/H^.., where /? is the approximate inverse 
phase transition duration and is the Hubble constant at the transition temperature. 
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nih = 35 GeV, = ^ 





FIG. 3: Calculated gauge dependence of phase transition parameters for a medium-mass Higgs 
boson. See fig. [2] for a thorough explanation of the different lines. 
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FIG. 4: Calculated gauge dependence of phase transition parameters for a high-mass Higgs boson. 
See fig. [2] for a thorough explanation of the different lines. 
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FIG. 5: Expected gravitational wave spectrum for a Higgs mass of 10 GeV, calculated in Landau 
gauge (.^ = 0), two high-^ gauges = 1,5), and the gauge-invariant Hamiltonian formalism. 
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FIG. 6: Expected gravitational wave spectrum for a Higgs mass of 35 GeV, calculated in Landau 
gauge = 0), one or two high-^ gauges = 1, 5), and the gauge- invariant Hamiltonian formalism. 
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FIG. 7: Expected gravitational wave spectrum for a Higgs mass of 120 GeV, calculated in Landau 
gauge (.^ = 0), a high-^ gauge = 1), and the gauge-invariant Hamiltonian formalism. 
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FIG. 8: Comparison of gravitational wave spectra calculated without daisy resummation in Landau 
gauge, and with resummation (dashed lines) in Landau gauge and two other gauges. 
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FIG. 9: Comparison of gravitational wave spectra calculated without daisy resummation in Landau 
gauge, and with resummation (dashed lines) in Landau gauge and one or two other gauges. 
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FIG. 10: Comparison of gravitational wave spectra calculated without daisy resummation in 
Landau gauge, and with resummation (dashed lines) in Landau gauge and one other R^ gauge. 
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